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We reveal the existence of high-density Feshbach resonances in the collision between the ground 
and metastable states of '^^^Yb and coherently produce the associated Feshbach molecules by pho¬ 
toassociation. The extremely small transition rate is overcome by the enhanced Franck-Condon 
factor of the weakly bound Feshbach molecule, allowing us to observe Rabi oscillations with long 
decay time between an atom pair and a molecule in an optical lattice. We also perform the precision 
measurement of the binding energies, which characterizes the observed resonances. The ultra-narrow 
photoassociation will be a basis for practical implementation of optical Feshbach resonances. 
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Remarkable advances in study on ultracold atomic 
gases owe much to their high controllability over various 
experimental parameters. A tunability of the interac¬ 
tions between atoms with Feshbach resonances (FRs) [1] 
is undoubtedly an important key to many experiments. 
One of the most significant applications of FRs is to 
create ultracold molecules. By adiabatically sweeping a 
magnetic field across a resonance, one can create ultra¬ 
cold, near-threshold Feshbach molecules, which plays a 
key role in many studies such as the crossover between a 
Bose-Einstein Condensate (BEC) and a Bardeen-Cooper- 
Schriefer superfluid [2-4], the initial step to form an 
ultracold, ro-vibronic ground state of (polar) molecules 
[5, 6], correlation measurements for ultracold atoms in an 
optical lattice [7-9], and so on. Photoassociation (PA) is 
another standard method to generate ultracold molecules 
which can create various molecules in electronic excited 
states with a simple optical excitation. [10]. However, 
the transition probability is never large enough even for 
a strong electric-dipole allowed (El) transition due to 
small overlap between the wave functions of free atoms 
and molecular bound states, resulting in small Franck- 
Condon factors. In addition, the radiative lifetime of 
the created molecular states is usually quite short for an 
El transition. As a consequence, coherent production of 
molecules by one-color PA has been impossible until re¬ 
cent demonstration with a BEC of ®®Sr [11], where the 
narrow ^5'o O ^Pi transition was exploited. 

In this Letter, we report on photoassociative creation 
of ultracold molecules associated with the ground ^S'o and 
long-lived metastable ^P 2 states of ^^^Yb with a radia¬ 
tive lifetime of longer than one second. Importantly, the 
otherwise quite small strength of the PA resonance is sig¬ 
nificantly enhanced by working around FRs for the closed 
channel. This enables us to demonstrate Rabi oscillations 
in the PA transition with a lifetime reaching hundreds 
of microseconds. Moreover, by using this narrow-line 
PA, observed FRs are characterized by precise measure¬ 
ment of the binding energies of near-threshold Fesbach 
molecules. Our observation of the strong PA line with 


a long expected radiative lifetime opens up the possibil¬ 
ity to suppress atom loss in optical Feshbach resonances 
[11-13]. 

First, we describe our observation of FRs in collisions 
between ^Sq and ^P 2 states of ^^^Yb, which is an impor¬ 
tant basis in this work. Most of the currently studied 
FRs are for alkali atoms, which originate from their hy- 
perfine structures with isotropic van der Waals interac¬ 
tions. On the other hand, the emergence of FRs associ¬ 
ated with the existence of anisotropic interactions was re¬ 
cently predicted [14] and observed in collisions of Er [15] 
and Dy [16]. For these atoms, anisotropy in an electro¬ 
static van der Waals potentials with non-zero electronic 
orbital angular momenta as well as dipole-dipole inter¬ 
actions induces couplings between an open channel and 
many closed channels with higher partial waves, lead¬ 
ing to an extremely high density of FRs [17]. In the 
present case of the collision between the ^Sq and ^P 2 
states of Yb atoms, anisotropy originates purely form 
the van der Waals interaction with an orbital angular 
momentum [18]. 

We start from the results of trap loss spectroscopy. 
Figure 1 (a) shows the relevant energy diagram and the 
experimental sequence. An ultracold Fermi gas of ^^^Yb 
with two spin components (jg-i-) = \^So{mi = -1-1/2)) 
and = \^So{mi = —1/2))) is prepared by sympa¬ 
thetic evaporative cooling with ^^^Yb in an optical dipole 
trap with a wavelength of 532 nm [19]. After evapo¬ 
ration we obtain a thermal gas of 8.3(2) x 10“^ atoms 
at a temperature T = 1.1(1) ^K. The trap frequen¬ 
cies are = (327,62.4,426) Hz with z in¬ 

dicating the direction of gravity, leading to the Fermi 
temperature of Tp = 620 uK. Remaining ^^^Yb atoms 
are removed by 556 nm light at the resonant frequency 
of the ^S'o —>■ ^Pi{F = 7/2) transition. At a static 
magnetic field of B = = 1.0 G, we excite a small 

fraction of atoms to the lowest hyperfine state in ^P 2 , 
|e) = \^P 2 {F = 3/2,wf = —3/2)) state, by applying 
resonant laser light at 507 nm [20]. Immediately af¬ 
ter the excitation, we ramp up the field to the desired 


2 


(a) 


(b) 


(c) 


(d) 




B[G] 


FIG. 1. (Color online), (a) Low-lying energy level diagram 
of (b-d) Observation of FRs (marked with rectangles) 

in a mixture of g^-g^-e, g^-e and respectively. Trap loss 
spectra at temperature of (b) 1.1 /rK, and (b,c) 1.8 /rK are 
shown. Pink points in (b) represents the data for a sample 
containing |e) atoms only. Red dashed line in (c) highlights 
the existence of a broad resonance at 6.6 G. 


value within 10 ms. After a hold time of 100 ms in the 
dipole trap, atoms in the ground state is removed, fol¬ 
lowed by repumping of survived atoms in the ^P 2 state 
into the ground state. The repumped atoms are detected 
by fluorescence imaging with a magneto-optical trap us¬ 
ing 399 nm cooling light. 

Figure 1 (b) shows 12 resonant losses which are at¬ 
tributed to the existence of FRs in \g'[)-\e) or |g 4 .)-|e) 
collisions [21]. We exclude the possibility of |e)-|e) reso¬ 
nances by repeating the same experiment with je) atoms 
only. Loss spectra for the |gt)-|e) and | 54 ,)-|e) states are 
also shown in Fig. 1 (c) and (d), respectively, by which 
we can identify the relevant spin combination for an open 
channel in each resonance. For spin polarization, we ap¬ 


ply optical pumping with the —>■ ^Pi(F = 1/2) tran¬ 

sition at the early stage of evaporative cooling. Due to 
the incompleteness of optical pumping, resonances with 
the residual spin component are also visible for each case 
of two-component spectrum. We estimate the ratio of 
atoms in each state to be (ffijUljs) = (30%, 50%, 20%), 
(70%, 10%, 20%) and (5%, 75%, 20%) for Fig. 1 (b), (c) 
and (d), respectively, which accounts for the contamina¬ 
tion of weak resonances. 

We note that the strong loss feature at 6.6 G seen in 
Fig. 1 (b) is the consequence of two overlapping reso¬ 
nances: the broad resonance in the jf/t)"!®) collision and 
the narrow resonance in the |<?^)-|e) collision. As shown 
in Fig. 1 (c), the former FR is accompanied by relatively 
small loss compared to the other narrow ones, which pos¬ 
sibly reflects the open-channel dominated feature of this 
resonance. The broad resonance with small loss is fa¬ 
vorable for practical use. Exploring FRs for fermionic 
isotopes is of great importance because the stability of 
resonantly interacting two-component Fermi gases pro¬ 
vides much wider practical applications. Especially in 
the case of resonances, using polarized fermions 

is the best way to avoid inelastic collisions in the ^P 2 
state which is a dominating decay channel [22]. 

By working around the observed FRs between the 
ground and metastable states of ^^^Yb, we directly form 
ultracold molecules with a photo-associative method with 
significantly enhanced rate. Owing to the ultranarrow 
O ^P 2 line in our PA method, it is possible to re¬ 
solve molecular resonance down to ~ 10 kHz, which is 
Doppler-limited linewidth of our samples. Figure 2 (a) 
shows a jgj^) je) laser spectrum taken at B — 6.2 G, 
in the vicinity of broad |5't-)"|6) resonance at 6.6 G. In 
this measurement, we simply irradiate excitation laser 
onto a ground state |gt)“l5i) mixture at T = 65(6) nK 
(T/Tp = 0.27) and measure the number of atoms remain¬ 
ing in the ground state by absorption imaging. On the 
red-detuned side of the atomic resonance, we clearly ob¬ 
serve a resonant loss of atoms, which can be attributed to 
the association of Feshbach molecules. Large spatial ex¬ 
tent of the Feshbach molecule provides the large Franck- 
Condon factor, giving PA rate comparable to the atomic 
transition. Later we return to this point in the context 
of atom-molecule Rabi oscillations in an optical lattice. 

The frequency detuning of the PA resonance from the 
atomic resonance directly gives the binding energy of 
the associated molecules. We systematically measure the 
magnetic field dependence of the binding energy and con¬ 
firm that the observed molecular state is associated with 
the FR at 6.6 G, as shown in Fig. 2 (b). In addition, we 
find another molecular level stems from the resonance at 
6.0 G and its anti level-crossing with one from the 6.6 G 
resonance (Fig. 2 (c)). For further characterization of 
the observed FR, we perform coupled two-channel cal¬ 
culation [23] of the binding energy. In the calculation, 
we fix the unknown parameters Cq = 3000 a.u. for the 























































3 


(a) 




5.6 6.0 6.4 

B[G] 



FIG. 2. (Color online), (a) ^So laser spectroscopy 

at B = 6.2 G. The PA resonance to the Feshbach molecule 
associated with the 6.6 G broadest FR is visible at —38 kHz 
from the atomic resonance. The PA laser pulse has a duration 
20 ms and intensity 2.3 W/cm^. (b) Binding energy of the 
Feshbach molecule associated with resonances at B = 6.0 G 
and 6.6 G. The solid lines are the result of the coupled two- 
channel calculation. Two dashed lines indicate the bare closed 
channel molecule and the shallowest bound state of the open 
channel. The dot-dashed line is the prediction of the universal 
relation Bb = —Ti^/ma^. (c) Detailed behavior near the anti 
level-crossing at 6 G. 


^So-^P 2 inter-atomic potential and the differential mag¬ 
netic moment between the open and the closed channel 
(5/i = 0.9/rB where /iB is the Bohr magneton. We find 
that the calculation with background scattering length 
Ubg = 290ao and the resonant width A = 1.11 G well re¬ 
produce the experimentally determined binding energy. 
In the low field B < 3 G, the data shows deviation 
from the theory, which indicates the existence of the cou¬ 
pling to unknown molecular levels. As seen in the figure, 
the binding energy of the Feshbach molecule remarkably 
deviates from that of the bare closed-channel molecule. 
Therefore the binding energy is not sensitive to the mag¬ 
nitude of Sfi, which prevents from extracting the mag¬ 
netic moment of the closed channel. 

To directly observe the induced change of the scatter¬ 
ing length, we perform laser spectroscopy on the ^S'o O 
^P 2 transition also in an optical lattice. In a deep op¬ 
tical lattice where each lattice site can be regarded as 
an isolated potential well, resonance from doubly occu¬ 
pied sites is shifted by the on-site interaction Ug^e — t^g,g 
which can be converted to the scattering length Qg^e 



FIG. 3. (Golor online), (a) Spectrum on the ^So ^P 2 
transition in an optical lattice with the depth of 14.7 Bb. (b) 
Magnetic field dependence of the scattering length de¬ 
termined by spectroscopy in an optical lattice. 


[18, 24]. We adiabatically load a degenerate spin mixture 
of jg-i-) and at T/Tp = 0.17(2) into an simple cubic 
optical lattice with a wavelength A = 532 nm. After 
reaching the lattice depth of 14.7 times the recoil energy 
Br = 7i^(27r/A)^/2m, we apply an excitation pulse with 
intensity 0.12 W/cm^ and duration 1 ms to excite \gi) 
atoms, followed by the repumping scheme to count |e) 
atoms. Due to the Pauli principle, only the combination 
of |(?t)"|6) is produced from doubly occupied sites and the 
corresponding scattering length ag^^e can be measured. 

Figure 3 (a) is a typical spectrum in the lattice, which 
shows three prominent peaks. The largest peak corre¬ 
sponds to the excitation of isolated atoms, accompanied 
by the blue sideband resonance separated by the band 
gap ^ 29 kHz. One more resonance found between these 
two resonances is the signal from the doubly occupied 
sites. We determine the scattering length ag^^e by using 
the analytic model of a harmonically trapped atom pair 
[25] as 


Oho r {-Ug„e/2nu^ - 1/2) 

^ r {-Ug^,e/2hu) ’ 

where the onsite interaction Ug^^e is given by the Ug^^e = 
2TThAf -I- Ug^^g^ with A/ denoting the frequency differ¬ 
ence between the singly and the doubly occupied sites 
and Ug^^g^^ = —2TTh x 82 Hz the onsite interaction in the 
ground state obtained by the known scattering length 
“ —0.15 nm [26]. The mean trap frequency 
uj/2tt = 28.3 kHz of each potential well is deduced from 
the band calculation [21]. 

We map out the magnetic field dependence of 
around the 6.6 G FR, as shown in Fig. 3. In addition 
to the resonant change of the scattering length at 6.6 G, 
the effect of the narrower 6.0 G resonance is also visi¬ 
ble. We fit the obtained field dependence with the model 
proposed in Ref [27], 


a{B) = Obg 
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FIG. 4. (Color online), (a) Rabi oscillations of isolated atoms 
(top) and atom-pairs molecules at various magnetic field 
(bottom) in an optical lattice. The laser intensity is fixed to 
58 W/cm^. For the atomic transition, the Rabi frequency 
is 11.8 kHz. (b) The atom-molecule Rabi frequency H as 
a function of the magnetic field. Symbols (circles, squares) 
distinguish the different branches of the anti-crossing at 6 G. 
The solid line shows the prediction from the wave function 
overlap (see text). 


where a = [(5^°^ — B^'^)j— B^'^ — 6 B 2 )]^. Here 
Bi) and B^^ are the zero-crossing points of bare 
molecular energies. For the shift SB 2 = B 2 — B^\ we 
adopt an approximation [23] 


SB2 — A2 


(Qbg/a)(l - abg/a) 
1 -I- (1 - Obg/a)^ 


(3) 


with the mean scattering length d = 4.4 nm evaluated 
from the Cq constant mentioned above. Equation 2 is 
valid when the width of the first resonance is much nar¬ 
rower than that of the second resonance, |Ai| <C IA 2 I. 
The fit yields Bi = 6.003(5) G, Ai = 3.7(3) mG, B 2 = 
6.607(9)=G, A 2 = 680(40) mG, and Obg = 406(15)ao. 
Obtained parameters slightly differ from that inferred 
from the binding energy measurement. Possible source 
of the uncertainty in the scattering parameters is the ex¬ 
istence of unknown resonances at higher fields [21]. Es¬ 
pecially, extracting the value of Obg from the fit is highly 
sensitive to tails of unknown resonances. 

Narrow linewidth of the PA transition and resulting 
long lifetime of an associated molecule enables coherent 
production of molecules. For the narrow line (F /2tt = 
7 kHz) PA of a ®®Sr BEG, coherent Rabi oscillations up to 
~ 10 fis were reported [11]. Isolated atom pairs residing 
on each site of an optical lattice are more suitable system 
to observe coherent phenomena, for unwanted inelastic 
collisions are suppressed and long lifetime of molecules 
can be achieved [28]. 

We demonstrate the Rabi oscillations between jg-t-)- 
Igj.) atom pairs and |5-t-)"|c) Eeshbach molecules via one- 
color PA. Figure 4 (a) shows both the atomic and atom- 
molecule Rabi oscillations in the optical lattice with 


lA-lEn depth. For the \gi) —>■ |e) transition is El al¬ 
lowed, the Rabi frequency is ^ 10 times larger than the 
purely M2 allowed case [29], which makes it much easier 
to observe clear oscillations. The e~^ decay time of the 
oscillations is on the order of 500 /iS, which is limited by 
the inhomogeneous broadening of the resonant frequency. 
The observed atom-molecule Rabi frequency is on the 
same order as that of the atomic transition and keeps 
reasonable magnitude over the range of several gauss. To 
account for the field dependence of the Rabi frequency, 
we carry out the following simple analysis. The optical 
coupling between the initial and the final state is propor¬ 
tional to {g^,e\d{R)\gf,gi), where d{R) is the transition 
dipole depending on the inter-atomic distance R [30]. For 
a weakly bound Feshbach molecule, d{R) is almost un¬ 
changed from its asymptotic value, i.e., the atomic tran¬ 
sition dipole, over the most range of molecular extent. 
Therefore the Rabi frequency can be reduced to the sim¬ 
ple product H = Hatom X \{g-f,e\g^,gi)\. Using the an¬ 
alytic expression for the two-particle wave functions in 
a harmonic potential [25] and the experimentally deter¬ 
mined scattering length, we find good agreement between 
the data and the above estimation (Fig. 4 (b)). 

In conclusion, we observed the FRs between the ground 
and metastable excited states of ^^^Yb and successfully 
formed associated Feshbach molecules by one-color PA. 
Observed density of resonances, 0.7 /G per spin combi¬ 
nation, is lower than the cases of highly magnetic Er and 
Dy, while much higher than the typical of alkali atoms. 
This may be the consequence of the simpler level struc¬ 
ture and less anisotropic interactions of ^^^Yb. Moder¬ 
ately dense resonances are convenient for tuning interac¬ 
tions with a selected well-behaved resonance and for fu¬ 
ture theoretical analysis of the observed resonances. The 
Feshbach molecule can be strongly coupled to the ground- 
state scattering wave function via the PA transition. Op¬ 
tical Feshbach resonance on this transition will effectively 
change the interaction between the ground state with a 
small atomic loss. Finally, we note that ultra-narrow PA 
demonstrated here is applicable not only to the homonu- 
clear collisions of alkaline-earth-like atoms, but also to 
the heteronuclear collisions between alkali atoms and 
alkaline-earth-like atoms, such as Li-Yb [31]. 
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SUPPLEMENTAL MATERIAL 

In this material, we use the following notation defined 
in the main paper: 

l5t) = =+1/2), (4) 

l5i) = =-1/2), (5) 

|e) = |3p2,F = 3/2,mF = -3/2). (6) 

EXCITATION TO THE STATE 

While both the Ig-j-) —>■ |e) and the Igj.) —>■ |e) transi¬ 
tions are allowed by the magnetic quadrupole (M2) se¬ 
lection rule, the latter transition is much stronger be¬ 
cause the hypefine-interaction induced electric dipole 
(El) transition is allowed [1]. In the low magnetic field, 
these transitions overlap each other and we chose the po¬ 
larization of the excitation laser to suppress the unwanted 
transition. However, when we excite jg-i-) atoms to the |e) 
state, the imperfection of the polarization creates a small 
number of jgj^) atoms by stimulated emission from the |e) 
state (see Fig. 5). In the experiment with polarized \gf) 
atoms (Fig. 1 (b) in the main paper), we apply the second 
pumping pulse after excitation to the |e) state, to remove 
Igj.) atoms. This operation heats up atoms by ~ 500 nK 
and is applicable only to non-degenerate thermal atoms. 
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FIG. 6. Trap loss spectrum for two-component mixtures of (a) 
|gi-t')-|e) and (b) \gi)-\e). (c) The same experiment with purely 
|e) atoms shows the absence of FRs within the |e) state. 


LOSS SPECTRUM AT HIGHER MAGNETIC 
FIELD 

In this section, we present the result of atomic loss 
spectroscopy of Feshbach resonances taken in an ex¬ 
tended field range. As in the case of Fig. 1 (b)-(d) in 
the main paper, we excite atoms to the |e) state, fol¬ 
lowed by magnetic-field sweep to variable final values at 
which atom loss is measured. For a technical reason, 
field ramp time after the state preparation is limited to 


30 ms in this range, which causes additional loss during 
the sweep. Figure 6 shows the 2 and 4 resonances for 
the |5t)"|6) l5i)“|c) mixtures, respectively. The res¬ 

onance at 11.5 G for the \g\)-\e) mixture shows broad 
loss which is comparable to the 6.6 G resonance studied 
in the main paper. Similarly, a broad resonance is also 
found for the \gi)-\e) mixture at 8.9 G. 

LIGHT SHIFT OF THE ^+2 STATE 
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Our 532 nm dipole trap and optical lattices provide the 
|e) state with slightly deeper potential than the \g) state. 
While the ac polarizability in the \g) state is independent 
of the polarization of the far off-resonant laser light, the 
light shift of the |e) state depends on the angle 9 between 
the polarization and the quantization axis [2, 3]. For 
linearly polarized laser light, the ac polarizability of a 
magnetic sublevel ttif in a hyperfine manifold F is given 
in terms of a scalar and a tensor polarizabilities as 


FIG. 5. (a) Optical Stern-Gerlach experiment for ^^^Yb po¬ 
larized in the Ig-i-) state by optical pumping before applying 
excitation laser (top) and immediately after excitation to the 
|e) state (bottom). Population in the \gi) state increases by 
the stimulated emission from the excited state, (b) Atom 
densities integrated along the vertical direction. 


= d"'’ + (3cos»9 - 1) 


(7) 

as long as the light shift can be treated as a small per¬ 
turbation to the Zeeman splittings [4]. 

Figure 7 shows the calculated ratio of the light shift 
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FIG. 7. The ratio of the light shift in the |e) state to that in 
the \g) state at 532 nm, as a function of the angle between 
the light polarization and the magnetic field. Here, we assume 
/og — 1.21 and = 0.17 for the scalar and tensor 

polarizability in the |e) state, respectively. 


in |e) state to that in the \g) state. In our setup, only a 
lattice beam along the x-axis has the polarization parallel 
to the magnetic field (z), giving rise to the maximum light 
shift for the ^P 2 state. Polarization of all other trapping 
beams, including the lattice beam along excitation laser 
(y), point to directions perpendicular to the magnetic 


field so that the differential light shift between the states 
\g) and |e) can be minimized. 

In order to determine the scattering length from spec¬ 
troscopy in an optical lattices (Fig. 3 of the main pa¬ 
per), we take account of the above light shift and the 
trap frequency lu of each potential well is calculated 
as y^^(w2'+'w 2)/2, where ojg = 2 tt x 26.6 kHz for the 
ground state is obtained by the band calculation and 
We = 27rx 29.9 kHz is the geometric average of the slightly 
anisotropic trap frequencies for the |e) state. Broaden¬ 
ing few kHz) of the linewidth of each peak is also 
caused by the differential light shift. In this respect, the 
instability of the magnetic field strength 0.5 mG), 
and the linewidth of the excitation laser (< 300 Hz) also 
contribute to the measured linewidth. 
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